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ABSTRACT 

Millimeter very-long baseline interferometry (mm- VLBI) provides the novel capacity to probe the emission 
region of a handful of supermassive black holes on sub-horizon scales. For Sagittarius A* (Sgr A*), the super- 
massive black hole at the center of the Milky Way, this provides access to the region in the immediate vicinity of 
the horizon. Broderick et al. (2009) have already shown that by leveraging spectral and polarization information 
as well as accretion theory, it is possible to extract accretion-model parameters (including black hole spin) from 
mm- VLBI experiments containing only a handful of telescopes. Here we repeat this analysis with the most recent 
mm- VLBI data, considering a class of aligned, radiatively inefficient accretion flow (RIAF) models. We find that 
the combined data set rules out symmetric models for Sgr A*'s flux distribution at the 3.9a level, strongly favoring 
length-to- width ratios of roughly 2.4: 1 . More importantly, we find that physically motivated accretion flow models 
provide a significantly better fit to the mm- VLBI observations than phenomenological models, at the 2.9cr level. 
This implies that not only is mm- VLBI presently capable of distinguishing between potential physical models for 
Sgr A*'s emission, but further that it is sensitive to the strong gravitational lensing associated with the propagation 
of photons near the black hole. Based upon this analysis we find that the most probable magnitude, viewing angle, 

and position angle for the black hole spin are a = 0.0 +0 64+0 - 86 , = 68°!2oo!2 8 ° > and f = - 52 °-i5°-24° east of north > 
where the errors quoted are the la and 2a uncertainties. 

Subject headings: black hole physics — Galaxy: center — techniques: interferometric — submillimeter — accre- 
tion, accretion disks 



1. INTRODUCTION 

Despite being invoked to power a variety of energetic as- 
trophysical phenomena, the detailed structure and dynamics 
of black hole accretion flows remain a central problem in as- 
trophysics. Moreover, using electromagnetic observations to 
probe the structure and dynamics of the black hole spacetimes 
requires a substantial understanding of the physical processes 
that determine the fate of the accreting matter. Only recently 
has it become possible to probe this physics via large-scale 
computational simulations. Nevertheless, ab initio calculations 
are beyond our present capability, requiring numerous simpli- 
fying, and in some cases unphysical, assumptions. This is evi- 
denced by the number of models proffered to explain the vari- 
ous properties of accreting black hole candidates. In turn, this 
ambiguity complicates efforts to use electromagnetic observa- 
tions to probe the structure and dynamics of the spacetime sur- 
rounding the black hole. 

By virtue of its proximity, the supermassive black hole at 
the center of the Milky Way, associated with the bright radio 
point source Sagittarius A* (Sgr A*), provides an unparalleled 
opportunity to study black hole accretion in detail. For this 
reason, Sgr A* may serve as an exemplar of the larger class 
of supermassive black holes specifically, and of black holes in 
general. Presently, the best estimates of the mass and distance 
of Sgr A* come from the observations of orbiting stars. These 
have yielded M = 4.3 ± 0.5 x 10 6 M© and D = 8.3 ± 0.4kpc, 
respectively, where both include the systematic uncertainties 
(Ghez et al. 2008; Gillessen et al. 2009b,a). The mass is nec- 
essarily confined to within the periapse of nearby stars, giving 
a maximum radius of roughly 10 2 AU ~3x 10 3 GM/c 2 , rul- 
ing out many extended objects. These represent the best mass 



measurement for any known black hole to date. 

In addition to the dynamical observations, a wealth of spec- 
tral and polarization data exists for Sgr A*. From these it is 
apparent that Sgr A* is unlike many active galactic nuclei, be- 
ing vastly underluminous, emitting a bolometric luminosity of 
roughly 10 36 erg, approximately 10~ 9 of Eddington. This is es- 
pecially small in light of the considerable amount of gas within 
the black hole's sphere of influence, presumably available for 
accretion (Loeb & Waxman 2007; Cuadra et al. 2008). As a re- 
sult it is widely accepted that Sgr A*'s accretion flow is qualita- 
tively different from those in its active analogs, though perhaps 
indicative of the roughly 90% of black holes that are presently 
not in an active phase. 

Nevertheless, the existing spectral and polarization data has 
produced a canonical set of components all models for Sgr 
A* include: populations of thermal and nonthermal electrons, 
nearly equipartition magnetic fields. Less certain is the struc- 
ture of the emission region. This is evidenced by the variety 
of models that have been proposed (e.g., Narayan et al. 1998; 
Blandford & Begelman 1999; Falcke & Markoff 2000; Yuan 
et al. 2002, 2003; Loeb & Waxman 2007). Despite being able 
to reproduce the observed features of Sgr A*, these differ dra- 
matically in the morphology of the emitting region. As a con- 
sequence, many of the theoretical ambiguities can be imme- 
diately addressed by direct probes of the spatial distribution 
of the emitting plasma surrounding the central supermassive 
black hole. 

The spectrum of Sgr A* peaks near millimeter wavelengths, 
implying a transition from optically thick to optically thin 
emission. The location of this emission is currently debated, 
however the presence of short-timescale variability at millime- 
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ter, near-infrared and X-ray wavelengths implies that optically 
thin emission is dominated by contributions arising in the im- 
mediate vicinity of the black hole. Furthermore, at millimeter 
wavelengths the blurring due to interstellar electron scattering 
is subdominant. Thus, at wavelengths of 1.3 mm and below it 
is possible to image the emitting region surrounding Sgr A*. 

Even with the strong gravitational lensing in the vicinity of 
the horizon, imaging the immediate vicinity of the black hole 
requires extraordinary resolutions. The silhouette cast by the 
horizon on the surrounding emission is roughly 53±2/ias 1 . 
At the present time, this resolution is accessible only via 
millimeter- wavelength very-long baseline interferometry (mm- 
VLBI). VLBI observations of Sgr A* at 1.4mm using the In- 
stitut de Radioastronomie Millimetrique (IRAM) 30 m tele- 
scope at Pico Veleta and one of the 15 m dishes at Plateau de 
Bure, produced the size estimate of 1 10 ± 60 //as, with the large 
uncertainties due to limited calibration accuracy (Krichbaum 
et al. 1998). 

The first successful mm- VLBI observation of Sgr A* with 
Earth-scale baselines was performed in April, 2007, during 
which visibilities were measured on the 4.6 x 10 3 km base- 
line between Mauna Kea, Hawaii to Mount Graham, Arizona 
(Doeleman et al. 2008). By fitting these with a gaussian model, 
Doeleman et al. (2008) found a typical intrinsic source size 
of 37^ /ias 2 (after correcting for the sub-dominant broadening 
due to interstellar electron scattering), smaller than the black 
hole silhouette. 

Since that time a number of groups have analyzed the 2007 
mm- VLBI data using various physically motivated accretion 
models for the emission region (Broderick et al. 2009; Huang 
et al. 2009; Moscibrodzka et al. 2009; Dexter et al. 2010), in- 
ferring from these efforts the black hole spin vector. Despite 
finding generally similar results, these have been limited by 
the lack of multiple long baseline observations and the limited 
north-south coverage obtained. Recently, a second, and con- 
siderably larger set of mm- VLBI observations have been re- 
ported (Fish et al. 2010), providing the opportunity to revisit, 
and substantially improve, constraints upon the black hole spin 
and accretion physics. 

Here we report upon the first effort to do this using a phys- 
ically motivated accretion model, similar to that described in 
Broderick et al. (2009), that fits the known spectral and polar- 
ization properties of Sgr A*. In addition to improving the re- 
sulting parameter estimation, it is now possible to identify sta- 
tistical signatures of both the asymmetry of the image and the 
importance of the underlying physics that governs the image 
morphology. Section 2 summarizes the full set of mm- VLBI 
observations we consider. Section 3 describes the models we 
consider and how the resulting visibility data is produced. How 
models are compared and the parameter estimates are produced 
is discussed in Section 4. The fitting process and results are 
presented in Section 5, and our best estimates for the black hole 
spin vector can be found in Section 6. Section 7 describes the 
implications for different potential future observations. Finally 
concluding remarks are collected in Section 8. 



1 The mass and distance measurements are strongly correlated, with mass 
scaling roughly as M oc D 1 8 (Ghez et al. 2008). 

2 We quote the la errors implied by the 3 a errors reported in Doeleman 
et al. (2008). 



2. SUMMARY OF MILLIMETER- VLB I OBSERVATIONS 

In the analysis presented here we make full use of the re- 
cent observations described in Fish et al. (2010) and Doeleman 
et al. (2008). In both cases, observations targeting Sgr A* were 
made at 1.3 mm using the Submillimeter Telescope (SMT) on 
Mt. Graham in Arizona, 10 m dishes in the Combined Ar- 
ray for Research in Millimeter- wave Astronomy (CARMA) at 
Cedar Flat, California, and the James Clerk Maxwell Telescope 
(JCMT) located on Mauna Kea, Hawaii. 

2.1. April 2007 

Doeleman et al. (2008) report upon measurements obtained 
on the nights of the April, 11 & 12, 2007, using the JCMT, 
SMT and a single CARMA dish. 19 visibility amplitudes 
were obtained on the CARMA-SMT and JCMT-SMT base- 
lines, with an upper limit on April 1 1th, 2007 along the JCMT- 
CARMA baseline. The locations of these observations on the 
u-v plane are indicated in the lower-left panel of Figure 1 , la- 
beled 2007. Signal-to-noise ratios typical of the short and long 
baselines are 8 and 4, respectively. 

During this time, observations the single-dish flux was esti- 
mated via the full CARMA array, operating as a stand-alone 
instrument, to be 2.4 ± 0.25 Jy. This is similar to the visibility 
amplitudes obtained on the CARMA-SMT baselines and con- 
sistent with a single, compact gaussian component (Doeleman 
et al. 2008). This flux is anomalously low in comparison to the 
typical 1.3 mm flux of ~ 3Jy, and was taken as evidence for 
Sgr A* appearing in a quiescent state. This interpretation is 
supported by the lack of a significant difference between anal- 
yses of each day separately (Broderick et al. 2009). 

Full details of the observations, calibration and data process- 
ing can be found in Doeleman et al. (2008). 

2.2. April 2009 

Fish et al. (2010) report upon more recent observations per- 
formed on the nights of April, 5-7, 2009, corresponding to the 
95, 96, and 97 days of 2009. These made use of the JCMT, 
SMT, and two CARMA dishes, operated as independent VLBI 
stations. 54 visibility amplitudes were obtained on JCMT- 
SMT and CARMA-SMT baselines on all days, and to both of 
the JCMT-CARMA baselines on days 96 and 97. Positions of 
the observations on each day are indicated in the upper panels 
of Figure 1, labeled 2009.95, 2009.96, and 2009.97. Signal-to- 
noise ratios typical of the short and long baselines are 17 and 
5, respectively. Thus, this second data set represents a signifi- 
cant improvement in both the number and precision of the data 
obtained. 

In addition to the VLBI baselines, the presence of two 
independent CARMA dishes in the array allowed the mea- 
surement of very- short baseline visibilities, probing angular 
scales ~ 10". These found substantially more correlated flux 
density than the CARMA-SMT baselines did, inconsistent 
with a single compact gaussian component. The interpreta- 
tion of the difference in correlated flux density between the 
inter-CARMA baselines and the CARMA-SMT baselines is 
presently unclear, and it may be possible for multiple geomet- 
ric models (e.g., annular rings, extended double source) to fit 
the data. Within the context of our analysis, we will assume 
that this difference is due to a separate large-scale component 
not present during the 2007 observations. This is indirectly 
supported by the fact that the source sizes inferred from the 
mid and long baseline data are unchanged despite the variations 
in the visibility magnitudes (Fish et al. 2010). Therefore, we 
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Figure 1. Locations in the w-v plane of the visibilities observed during the 2007, 2009.95, 2009.96, 2009.97 epochs. Also shown are the combined set of observations. 
Finally, for references the combined set is compared to the potential baselines from existing and upcoming sub-mm telescopes. Each baseline is color-coded according 
to the associated two sites. In all plots, detections are denoted by green circles and upper-limits by red triangles. 



do not consider the inter-CARMA data further here, restricting 
ourselves to modeling the compact component observed with 
the longer baselines. 

On days 95 and 96 the short-baseline flux densities are con- 
sistent with each other, with inferred single dish fluxes of 
2.15 ±0.06Jy, which while somewhat lower than those ob- 
tained in 2007, justify treating these as a similar quiescent 
period. This is not the case for day 97, which exhibited a 
30%-40% increase in the luminosity of the compact compo- 
nent. Note that during the 3 hr observing periods on days 96 
and 97 there is no evidence for rapid changes in the CARMA- 
SMT visibility amplitudes, implying that during each Sgr A* 
was stable; i.e., the process responsible for the brightening oc- 
curred between observing periods and is stable on timescales 
of hours. As a consequence, we will treat the visibilities ob- 
tained on each day as due to a stationary source, though with 
properties that vary from day to day. 

Full details of the observations, calibration and data process- 
ing can be found in Fish et al. (2010). 

2.3. Combined Data Set 

Combined, the 2007 and 2009 mm- VLB I measurements may 
be separated into 4 observational epochs: that containing the 
entire 2007 observations (2007), those on day 95 of 2009 
(2009.95), those on day 96 of 2009 (2009.96), and those on day 
97 of 2009 (2009.97). The combined coverage in the u-v plane 
is shown in the lower-middle panel of Figure 1 . The long base- 



lines (JCMT-CARMA and JCMT-SMT) are oriented primar- 
ily in the east- west direction, extended roughly 3.6 GA. Never- 
theless, the combined data set also extends roughly 2GA in the 
north-south direction, providing substantial angular coverage 
in the u-v plane for the first time. 

In Section 7 we will discuss the implications out analysis 
has for future observations. However, we note here that the 
baselines considered in the 2007 and 2009 mm- VLB I exper- 
iments are a small subset of the baselines that are possible 
with existing mm and sub-mm telescopes. Figure 1 shows 
the combined visibility data set in comparison to baselines as- 
sociated with other potential mm- VLB I stations. These in- 
clude stations in Chile (e.g., the Atacama Pathfinder Experi- 
ment, Atacama Submillimeter Telescope Experiment, and At- 
acama Large Millimeter Array; APEX, ASTE, and ALMA, re- 
spectively), Mexico (Large Millimeter Telescope; LMT), Spain 
(Pico Veleta; PV), France (Plateau de Bure; PdB), and at 
the South Pole (South Pole Telescope; SPT). These both, ex- 
tend the region covered in the u-v plane, and provide addi- 
tional complementary short and intermediate baselines, primar- 
ily along the north- south directions. To date, visibilities on 
only a handful of potential baselines have been measured. 

3. VISIBILITY MODELING 

Our primary goal is to use physically motivated models of 
Sgr A*'s accretion flow to infer the properties of the central 
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supermassive black hole and its surrounding matter. To do this 
we compare both physical and phenomenological models of 
Sgr A* to the mm-VLBI visibilities. This requires the compu- 
tation of model visibilities. Given a trial image intensity dis- 
tribution, /(a,/?), where a and f3 are angular coordinates, we 
may compute the visibilities in the standard fashion: 



V(u,v) = 



dadPe- 27ri(au+(3v)/x I(a,f3). 



(1) 



Here we describe three classes of model images: those asso- 
ciated with radiatively inefficient accretion flows (RIAFs) of 
the form discussed in Broderick & Loeb (2006a), symmetric 
and asymmetric gaussians. We also summarize the effects of 
interstellar electron scattering. 

3.1. Radiatively Inefficient Accretion Flows 

We employ a suite of radiatively inefficient accretion flow 
(RIAF) models, first described in Broderick & Loeb (2006a), 
and based upon those of Yuan et al. (2003). Here these models, 
which henceforth we refer to as BL06, are summarized. 

Sgr A* transitions from an inverted, presumably optically 
thick spectrum to an optically thin spectrum near millimeter 
wavelengths. This implies that near 1.3mm Sgr A* is only be- 
coming optically thin, and thus absorption in the surrounding 
medium is likely to be important. This transition does not oc- 
cur isotropically, happening at longer wavelengths for gas that 
is receding and at shorter wavelengths for gas that is approach- 
ing. Therefore, properly modeling the structure and relativistic 
radiative transfer is crucial to producing high fidelity images. 

Although Sgr A* is vastly sub-Eddington, its bolometric lu- 
minosity, roughly 10 36 ergs _1 , is still large in absolute terms, 
Like many AGN, in the radio Sgr A* exhibits the nearly-flat, 
power-law spectrum associated with non-thermal synchrotron 
sources, with the power emitted (vL u ) peaking at millimeter 
wavelengths. As a consequence, it has been widely accepted 
that Sgr A* is accretion powered, implying a minimum ac- 
cretion rate of lO~ lo M0yr _1 . It is presently unclear how this 
emission is produced, evidenced by the variety of models that 
have been proposed (e.g., Narayan et al. 1998; Blandford & 
Begelman 1999; Falcke & Markoff 2000; Yuan et al. 2002, 
2003; Loeb & Waxman 2007). Models in which the emis- 
sion arises directly from the accreting gas have been subsumed 
into the general class of RIAFs, defined by the generally weak 
coupling between the electrons, which radiate rapidly, and the 
ions, which efficiently convert gravitational potential energy 
into heat (Narayan et al. 1998). This coupling may be suffi- 
ciently weak to allow accretion rates substantially in excess of 
that required to explain the observed luminosity with a canon- 
ical AGN radiative efficiency of 10%. However, the detec- 
tion of linear polarization in Sgr A* above 100GHz (Aitken 
et al. 2000; Bower et al. 2001, 2003; Marrone et al. 2006) 
and subsequent measurements of the Faraday rotation mea- 
sure (Macquart et al. 2006; Marrone et al. 2007), have implied 
that the accretion rate near the black hole is much less than 
the Bondi rate, requiring the existence of large-scale outflows 
(Agol 2000; Quataert & Gruzinov 2000). 

Relating the outflow to the properties of the accretion flow 
requires an ab initio calculation that is presently not possi- 
ble. Nevertheless, a number of authors have studied this re- 
lationship in the context of a variety of simplifying assump- 
tions, with large-scale general-relativistic magnetohydrody- 
namic and radiative-hydrodynamic simulations playing a cen- 
tral role (De Villiers et al. 2005; McKinney 2006; Hawley & 
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Figure 2. Comparison of the spectrum of the most probable accretion model 
and the observed SED of Sgr A*. Orange circles are from Yuan et al. (2004), 
and references therein, for which the errorbars are indicative of the variabil- 
ity. Yellow squares are coincident flux measurements from Marrone (2006), 
for which the errorbars are indicative of the intrinsic measurement error. The 
green bar shows the flux range of the compact component inferred from fitting 
accretion models to the mm-VLBI data (Section 5. In addition to the most 
probable model, spectra are shown for emission from a rapidly rotating black 
hole (a = 0.998) as seen nearly face on (0 = 1°, green dotted) and edge on 
(6 = 90°, purple dotted), indicating the range of variation within the image li- 
brary. Finally, for reference the contribution from the thermal (red dash) and 
nonthermal (blue long-dash) are shown. Note in particular that at 1.3mm the 
nonthermal contribution is not negligible. 

Krolik 2006; Beckwith et al. 2008; McKinney & Blandford 
2009; Tchekhovskoy et al. 2010; Dexter et al. 2010; Penna 
et al. 2010; Kurosawa & Proga 2009). In these it has been 
found that the structure and dynamics of the outflow critically 
depends upon the initial conditions. The applicability of the 
MHD prescription to Sgr A* is still unclear, where the accre- 
tion rate is sufficiently low that non-MHD effects may become 
important (Sharma et al. 2006, 2007). More importantly, most 
of these approaches do not model the electron heating (beyond 
ad hoc prescriptions) and none model the production of non- 
thermal electrons (see, e.g., Moscibrodzka et al. 2009; Dexter 
et al. 2010; Shcherbakov et al. 2010). Furthermore, simulations 
are computationally expensive to produce. For these reasons 
we adopt a simple, self-similar model for the accretion flow 
which includes substantial mass loss. 

For concreteness, as in Broderick & Loeb (2006a), we follow 
Yuan et al. (2003) and employ a model in which the accretion 
flow has a Keplerian velocity distribution, a population of ther- 
mal electrons with density and temperature 



rs 



-l.i 



~Z 2 /2P 2 



and 



rs 



-0.84 



(2) 



(3) 



respectively, and a toroidal magnetic field in approximate (J3 = 
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Figure 3. Best fit images and visibilities for the three models considered: symmetric gaussian (top), asymmetric gaussian (middle) and BL06 accretion flow (bottom). 
For each model we show the intrinsic flux distribution (left), flux distribution after interstellar electron scattering (center) and the visibility amplitudes (right). For 
reference the locations of the observed visibilities (over all epochs) are shown by the white points. In all plots the intensity scales linearly with / and V. 



10) equipartition with the ions (which are responsible for the 
majority of the pressure), i.e., 



8tt 



= P n e ,± — 7z — 
12r 



(4) 



In all of these, rs = 2GM/c 2 is the Schwarzschild radius, p is 
the cylindrical radius and z is the vertical coordinate. Inside of 
the innermost-stable circular orbit (ISCO) we assume the gas is 
plunging upon ballistic trajectories. In principle the plunging 
gas can still radiate, though in practice it contributes little to the 
overall emission due to the large radial velocities it develops. 
In the case of the thermal quantities the radial structure was 
taken from Yuan et al. (2003), and the vertical structure was 



determined by assuming that the disk height is comparable to 
p. Note that all of the models we employ necessarily have the 
spin aligned with the orbital angular momentum of the accre- 
tion flow. For the regions that dominate the mm emission, this 
assumption is well justified due to disk precession and viscous 
torques, though it may be violated at large distances. 

Thermal electrons alone are incapable of reproducing the 
nearly-flat spectrum of Sgr A* below 43 GHz. Thus it is nec- 
essary to also include a nonthermal component. As with the 
thermal components, we adopt a self- similar model for a pop- 
ulation of nonthermal electrons, 



rs 



-2.02 



~z 2 /2p 2 



(5) 
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with a power-law distribution corresponding to a spectral in- 
dex of 1 .25 and cut off below Lorentz factors of 10 2 (consistent 
with Yuan et al. 2003). The radial power-law index was chosen 
to reproduce the low frequency spectrum of Sgr A*, and is in- 
sensitive to the black hole properties due to the distant location 
of the long-wavelength emission. 

The primary emission mechanism at the wavelengths of in- 
terest is synchrotron, arising from both the thermal and non- 
thermal electrons. We model the emission from the thermal 
electrons using the emissivity described in Yuan et al. (2003), 
appropriately altered to account for relativistic effects Broder- 
ick & Blandford (see, e.g., 2004). Since we perform polar- 
ized radiative transfer via the entire complement of Stokes pa- 
rameters, we employ the polarization fraction for thermal syn- 
chrotron as derived in Petrosian & McTiernan (1983). In doing 
so, we have implicitly assumed that the emission due to ther- 
mal electrons is isotropic, which while generally not the case is 
unlikely to change our results significantly. For the nonthermal 
electrons, we follow Jones & O'Dell (1977) for a power-law 
electron distribution, with an additional spectral break associ- 
ated with the minimum electron Lorentz factor. For both emis- 
sion components the absorption coefficients are determined di- 
rectly via Kirchhoff 's law. Images are then produced using the 
fully relativistic ray-tracing and radiative transfer schemes de- 
scribed in Broderick & Loeb (2006a,b) and Broderick (2006). 
An example image and associated visibilities are shown in the 
bottom line of Figure 3. 

Because Yuan et al. (2003) neglected relativistic effects and 
assumed spherical symmetry, it is not directly applicable here. 
For these reasons, as in Broderick & Loeb (2006a), the coef- 
ficients (n® th ,T®,n® nth ) were adjusted to fit the radio spectral 
energy distribution (SED) of Sgr A*, shown in Figure 2. The 
1.3 mm and 0.43 mm points (orange squares) were measured 
coincidentally, and the errors represent the intrinsic measure- 
ment error (Marrone 2006). All other points (red circles) are 
taken from Yuan et al. (2004) (and references therein) and were 
obtained by averaging over multiple epochs. As a result the er- 
rorbars represent the range of variability, and are correspond- 
ingly larger. 

As in Broderick et al. (2009) we systematically fit 
Sgr A*'s SED at a large number of positions in the 
spin-inclination (a-0, where is the viewing angle rel- 
ative to the disk axis) parameter space, specifically at 
values of AG {0,0.1,0.2, ...,0.9,0.99,0.998} for all e 
{ 1° , 1 0° , 20° , . . . , 80° , 90° } , producing a tabulated set of the co- 
efficients (Hgth>^?°> n <?nth) at 120 points in the a-6 parameter 
space. In all cases it was possible to fit the SED with extraordi- 
nary precision, with reduced x 2 < 0.3 in all cases and typically 
reduced x 2 — 0.17. This is likely a consequence of employ- 
ing the variability-determined errorbars on the non-coincident 
flux measurements. Over the a-6 plane the x 2 was remarkably 
uniform, implying that on the basis of the spectra alone it is 
difficult to constrain the parameters of our simple model. Nev- 
ertheless, the dynamical properties of the disk manifest them- 
selves in the breadth of the sub-millimeter bump. Models with 
low spin and/or low 6 have little Doppler shifting, and corre- 
spondingly narrow bumps. In contrast, models with large spins 
(a > 0.9) viewed edge-on (6 > 80°) have broad bumps, and 
are responsible for the relatively larger, though still small, x 2 - 
From the tabulated values, the coefficients are then obtained at 
arbitrary a and 6 using high-order polynomial interpolation. 

An example spectrum, resulting from the above procedure, 
is shown in Figure 2, corresponding to a = 0.01 and 6 = 66°. 



In addition, the individual contributions from the thermal and 
nonthermal components are also shown (though the cross- 
absorption is neglected). The necessity of the nonthermal elec- 
trons is clearly illustrated at long wavelengths, where the ther- 
mal contribution is negligible. The transition from nonther- 
mally dominated to thermally dominated occurs near 2 mm, 
though the precise location depends upon a and 6. However, 
note that at no point can the nonthermal component be ne- 
glected. Specifically, at 1.3 mm the nonthermal component is 
still responsible for roughly 30% of the emission. As a re- 
sult, efforts to model the millimeter image of Sgr A* without 
accounting for the nonthermal component can produce order 
unity systematic errors in parameter estimation. 

During the mm-VLBI observations Sgr A*'s flux varied by 
roughly 30%. We model this as a variable accretion rate, 
moving the electron density normalization up and down. In 
practice, we reduced the electron density normalization by an 
amount sufficient to produce a total flux of 2.5 Jy, and then 
multiplied the resulting images by a correction faction dur- 
ing the mm-VLBI data analysis. Because the source is not 
uniformly optically thin, this is not strictly correct, though 
this makes a small change to the images themselves. For 
the purpose of the mm-VLBI data analysis (described below) 
we produced 9090 images, with flux normalized as described 
above, at a e {0,0.01,0.02, ...,0.98,0.99,0.998} for each 6 e 
{ 1 ° , 2° , . . . , 89° , 90° } . We then produce models with arbitrary 
position angles, £, by rotating the image on the sky. For this 
purpose we define £ such that at 0° the projected spin vector 
points north, and as £ increases points progressively more east- 
ward 3 . 

The Faraday rotation measures observed in Sgr A* are pro- 
duced at much larger radii than those of interest in direct imag- 
ing experiments. Nevertheless, it is worth noting that the mod- 
els employed here are broadly consistent with the polarization 
observations, though breaks in the radial power-laws which de- 
fine the properties of the thermal electron component may be 
required at large spins. 

For the purposes of fitting the mm-VLBI visibilities, during 
each observation epoch this model has 4 parameters: spin (a), 
viewing angle (6), position angle (£) and flux normalization 
(Voo). When we analyze multiple epochs together the parame- 
ters defining the orientation of the system (a,0,£) will be held 
fixed, while those corresponding to the time- variable accretion 
rate (Voo) will be allowed to vary, though a full discussion will 
have to await Section 5. 

3.2. Gaussian Flux Distributions 

For comparison we consider two gaussian flux distributions. 
These differ from the accretion flow model described above in 
that they are purely phenomenological, without any clear phys- 
ical motivation and thus not constrained at all by the spectral 
and polarization properties of Sgr A*. As a result, we might 
expect these to be intrinsically less likely than physically mo- 
tivated models that are already chosen to be consistent with 
these properties. Nevertheless, we will consider them on equal 
footing with the BL06 model discussed above. For reasons that 
will become clear, we consider both symmetric and asymmet- 
ric gaussian intrinsic flux distributions. 

We may write the asymmetric gaussian flux distribution as 

3 Note that this is opposite the definition employed in Broderick et al. 
(2009). 
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where ciM^m and OM,m are the angular coordinates and widths in 
the major/minor axis directions. This is fully defined once aM, m 
and the position angle of the major axis is given. However, we 
choose to parametrize the asymmetric gaussian in terms of a 
single width, a, an anisotropy parameter, A, and the position 
angle: 



a 2 a 2 
I = Voo exp ( - cos 205 



(7) 



where a = and 03 is the angular coordinate measured 

from the position angle. The a and A are related to the a M ,m by 



1 1 1 

a 2 2a 2 2al 



and 



1 



1 



M 



t 2 2a 2 m 2a 2 M ' 



(8) 



Example gaussian images, with associated visibilities, are 
shown in the top two lines of Figure 3. 

Clearly, for isotropic configurations (i.e., cjm = <J m ), A = and 
& = &M,m- More generally, cr m /cr M = y/(l -A)/(l +A). Thus, 
this parametrization has the virtue of separately emphasizing 
size (via a) and asymmetry (via A) in the image. Note that 
this model has precisely the number of free parameters as the 
accretion flow model described above: those describing the im- 
age morphology, (<r,A,£), and the flux normalization, Voo, for 
each epoch. 

While the symmetric gaussian models are obviously a subset 
of the asymmetric gaussian models (corresponding to when A 
vanishes), we must be careful to distinguish the number of free 
parameters. In this case, A and £ are superfluous, and for each 
epoch we have only 1 parameter. 

3.3. Interstellar Electron Scattering 

The effect of interstellar electron scattering in the direction 
of Sgr A* have been carefully characterized empirically by a 
number of authors. This has been found to be consistent with 
convolving the source with an asymmetric gaussian, with major 
axis nearly aligned with east- west, and a A 2 wavelength depen- 
dence. We employ the model from Bower et al. (2006), which 
has major axis oriented 78° east of north, with associated full 
width at half-maximum for the major and minor axes given by 



FWHM^ S = 1.309 



A 
lcm 



FWHM|f = 0.64 



lcm 



mas, 



mas, 



(9) 



respectively. In practice, the interstellar electron scattering 
convolution was effected in the u-v plane, where the convo- 
lution reduces to a multiplication. 

4. BAYESIAN DATA ANALYSIS 

In fitting the observed visibilities we have two primary goals: 
choosing among various possible model flux distributions and 
estimating the parameters of these models. Both of these may 
be naturally accomplished within the context of Bayesian anal- 
ysis. Here we briefly summarize how we do this. 

We define the likelihood, p(\\q), for observing the visibili- 
ties V given the model parameters q, as described in Broderick 
et al. (2009). From this we obtain the log-likelihood, which we 
refer to as x 2 '- 

X 2 =-21ogMV|q) + C, (10) 



where the normalization constant depends only upon the partic- 
ulars of the data and, since we will only be interested in com- 
paring identical data sets, will henceforth be ignored. When 
only detections are considered this reduces trivially to the stan- 
dard definition of x 2 . Here it differs only due to the pres- 
ence of an upper-limit upon the visibility along the CARMA- 
JCMT baseline during the 2007 epoch. 

4.1. Model Comparison 

To compare the significance of different models we make use 
of the Bayesian Information Criterion (BIC) and the Akaike In- 
formation Criterion (AIC). Both of these are discussed in de- 
tail within the context of astrophysical observations by Liddle 
(2007) and Takeuchi (2000), and references therein. Thus, here 
we only define and summarize the properties of these statistics. 

In terms of the smallest effective x 2 f° r a given model, the 
BIC is defined by 



BIC = xL + *ln#, 



(11) 



where k is the number of model parameters and N is the number 
of data points. Note that this is simply the x 2 statistic penal- 
ized for models with large numbers of parameters. Assuming 
that the data points are independent (likely true) and identically 
distributed (nearly true), this is related to the posterior proba- 
bilities of two different models, M\ and M 2 , by 



p(Mi\\) _ p(Mj) 
p(M 2 \Y) p{M 2 ) i 



-(BICi-BIC 2 )/2 



(12) 



where p(M\,i) is me prior on model M\^ 2 . Thus, up to the un- 
known priors, the BIC is a measure of the posterior probability 
for a given model. If we further assume that p(M\) = p(M 2 ), 
the difference in BIC's gives the relative posterior probabilities 
directly. Therefore the model with the lowest BIC is preferred. 
Similarly, the AIC is defined by 



AIC = xL + 2*+ 



2k(k+ 1) 
N-k-1 



(13) 



where we have included a correction appropriate for when N is 
small (Burnham & Anderson 2002, 2004). This is very simi- 
lar to the BIC: the x 2 statistic penalized by a factor depending 
upon the number of model parameters, though with a some- 
what different penalty. Unlike the BIC, the AIC is not directly 
related to the posterior probability of a given model. Rather it is 
an approximate measure of the difference between the true data 
distribution and the modeled data distribution. Nevertheless, it 
is possible to interpret the AIC in terms of a model likelihood 
in a way identical to the BIC. 

For either criterion, lower values are preferred, with the rel- 
ative significance given by 



w u =e 



-(IQ-IC;)/2 



(14) 



where IC may be replaced by BIC or AIC. The AIC are con- 
ventionally judged on the Jeffrey's scale, which sets AIC > 5 
as "strong" and AIC > 10 as "decisive" evidence against the 
model with the higher IC. However, here we describe these in 
terms of the typical a as well, with model i being excluded at 
the ncr level if there exists a model j for which w^ is less than 
the associated cumulative normal probability (e.g., la implies 
that Wij < 0.32, 2a implies that w^ < 0.05, 3 a implies that 
Wij < 0.003, etc.). 
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4.2. Parameter Estimation 

For the accretion flow model we have the additional problem 
of identifying the most likely model parameters. The procedure 
we use to estimate the posterior probabilities is identical to that 
described in Broderick et al. (2009). In particular, we assume 
flat priors on all of the visibility normalizations, Voo, the spin 
magnitude, a, and an isotropic prior on the spin direction (0, £). 
Since we are primarily interested in the estimates for the black 
hole spin, we present the posterior probability of a marginal- 
ized over the Voo, We also construct marginalized poste- 
rior probability distributions of a, 6, and £ in the normal way 
(for specific definitions of p(a), p(0), and /?(£)> see Broderick 
et al. 2009). 

5. MODEL FITTING 

A number of important implications follow from computa- 
tions of the relevant x 2 f° r me three image models described 
in Section 3. These include whether or not we are justified in 
comparing mm-VLBI observations obtained at different times, 
the symmetry of the image and the importance of physics for 
reproducing the measured visibilities. 

5.7. Consistency of the 2007 & 2009 Epochs 

The dynamical timescale of the Sgr A* is comparable to the 
orbital period at the ISCO, as measured at Earth, and ranges 
from roughly 4 min to 30 min, depending upon spin. As a re- 
sult it is not at all clear that we may ignore the possibility of 
variability when attempting to model mm-VLBI observations 
spanning many nights, let alone years. Broderick et al. (2009) 
took special pains to ensure that the visibilities measured on the 
two consecutive days were consistent with a single, static un- 
derlying flux distribution. This was done by comparing fits to 
the individual days. Here we repeat this analysis for the 2009 
observations as well. 

In 2007, at 2.4 Jy, the luminosity of Sgr A* was anomalously 
low and stable over the two observation days. This is not the 
case during the 2009 observation, during which Sgr A* exhib- 
ited a dramatic brightening on the third day. During the preced- 
ing two days the luminosity of Sgr A*'s compact component, 
corresponding to scales smaller than 10 2 rs, was significantly 
smaller than that associated with the 2007 observations. Thus 
it is clear from the outset that we are not justified in compar- 
ing a single, static model to the observations. Instead, we be- 
gin with the ansatz that the morphology of Sgr A*'s image is 
fixed, with the flux variability being driven by changes in the 
accretion rate on day-to-day timescales. While this period is 
considerably larger than the 30 min timescale over which the 
properties of the accretion flow may change, it is justified in 
part by the stability of Sgr A*'s luminosity on these scales as 
well as the intrinsically short duration (~ 2hr) of the observa- 
tions each night. Therefore, we separate the data into 4 epochs: 
2007, 2009.95, 2009.96 and 2009.97, corresponding to the data 
obtained in 2007 and on days 95, 96 and 97 of 2009, respec- 
tively. For all epochs we keep the parameters that define the 
image morphology fixed, e.g., (cr,A,£)» or but al- 

low the overall flux normalization to vary from epoch to epoch. 
Upon fitting each epoch separately, and all epochs together, we 
may ask if the resulting parameter likelihoods are consistent 
with each other, i.e., check if our ansatz is self consistent. We 
will remark upon this further in the sections describing the fits 
for the individual models; however here it is sufficient to note 
that in all cases we find that the epochs are indeed consistent 
with a single underlying image morphology. 
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Figure 4. Symmetric gaussian x 2 's as a function of a, at the most likely Voo, 
for the 2007, 2009.95, 2009.96, and 2009.97 epochs, as well as when all epochs 
are combined. The minimum % 2 , number of visibility observations (N) and fit 
parameters (k) are listed in each plot. In all cases the range shown corresponds 
to min(x 2 )- 1 to min(x 2 )+6. For references the la error estimate is shown by 
the dotted errorbars for each data set. 

5.2. Symmetric Gaussian 

We begin with the symmetric gaussian model, which is pri- 
marily sensitive to the characteristic size of the image. In this 
case for each epoch there are 2 parameters: a and Vbo- The 
minimum x 2 is shown as a function of a for each epoch in Fig- 
ure 4. With the exception of epoch 2009.96, the reduced-% 2 
is comparable to unity. Over all epochs the a with the highest 
likelihood varies over 1.5 /ias, well within the single-epoch un- 
certainty (defined by the region in which x 2 is within unity of 
the minimum value). In particular, there are no trends distin- 
guishing either the much earlier 2007 epoch or the considerably 
brighter 2009.97 epoch. Thus we conclude that the character- 
istic size of Sgr A* did not vary substantially from one epoch 
to the next, despite considerably changes in its luminosity. 

The best fit intrinsic source size is a = 15.8 ±0.2 //as 
(FWHM = 37.2 ± 0.5 /ias) with the flux normalization provided 
in Table 1, ranging from 2.07 Jy to 2.81 Jy over the various 
epochs. The associated reduced-% 2 is 1.16, with 66 degrees 
of freedom. Upon convolving with the interstellar electron 
scattering, which along the general direction of the baselines 
SMT-JCMT and CARMA-JCMT baselines has a FWHM of 
width of 22 /ias, we find the FWHM of the broadened image is 
43.2 dz 0.6 /ias. This is in excellent agreement with the inferred 
size of 43^ found by Doeleman et al. (2008) on the basis of 
the 2007 epoch alone 4 . It is also in excellent agreement with 
the inferred sizes of 41.3+}^, 44.4+}, and 42.6!} found by Fish 
et al. (2010) for epochs 2009.95, 2009.96, and 2009.97, respec- 
tively 4 . The associated intrinsic and scatter-broadened image, 
with the associated visibilities are shown in the left, center and 

4 Here we have quoted the la errors upon their result in order to provide a 
direct comparison. 
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Table 1 

Model Fitting Results Summary 



Model k x 2 DoF X 2 /DoF V 2 007 V 2 009 - 95 V 2 009 - 96 V™ 9 - 91 BIC wff w6 AIC w£ 



Estimated Errors 

Symmetric Gaussian 5 76.78 66 L16 237 Z08 Z03 Z88 90 5 x 1(T 4 87/7 8 x 1(T 5 

Asymmetric Gaussian 7 61.50 64 0.961 2.53 2.25 2.23 3.06 91.3 1 x 10~ 2 77.3 1 x 10~ 2 

BL06 7 53.09 64 0.830 2.45 2.18 2.16 3.00 82.9 1 68.9 1 



Implied Errors 

Symmetric Gaussian 5 92.51 66 L40 237 2M Z03 Z88 114 4 x 10~ 5 103 9 x 10~ 6 

Asymmetric Gaussian 7 74.10 64 1.16 2.53 2.25 2.23 3.06 104 6 x 10 -3 89.9 6 x 10 -3 

BL06 7 63.97 64 1.00 2.45 2.18 2.16 3.00 93.8 1 79.7 1 



right panels of the top row of Figure 3. 

5.3. Asymmetric Gaussian and Image Anisotropy 

During the 2007 epoch the mm- VLB I visibilities are con- 
centrated nearly exclusively along a single line, oriented nearly 
east- west (see the top-left panel of Figure 1). However, upon 
including the 2009 epochs, the portion of the u-v plane sam- 
pled covers a bow-tie shaped region with opening angle 26° 
(see the bottom-center panel of Figure 1). While this is insuffi- 
cient to generate an image directly, the coverage is sufficient to 
address the gross angular structure of Sgr A*'s image. Where 
the symmetric gaussian provides a phenomenological way in 
which to estimate the typical size of Sgr A*'s emitting region, 
an asymmetric gaussian can begin to probe its symmetry. 

Figure 5 shows the minimum x 2 ( or equivalently, the max- 
imum likelihood) as a function of the average size, a, and 
anisotropy parameter, A. The four left panels show this for the 
individual epochs, while the large right panel shows this for the 
combined data set. Unlike the symmetric gaussian model, the 
likely regions have somewhat different morphologies. This is 
due to the different coverage of the u-v plane during the var- 
ious epochs (for example, the likely regions for epochs 2007 
and 2009.95 are similar because the u-v coverage during these 
epochs is similar). Nevertheless, the region preferred by the 
combined data sets is present in all cases, implying that as with 
the symmetric gaussian all epochs are consistent with a sin- 
gle underlying image morphology. During this time the flux 
normalization of the compact component varied from 2.23 Jy 
(2009.96) to 3.06 Jy (2009.97). 

For all epochs the reduced-% 2 is nearly unity, ranging from 
0.45 (2007) to 1 .34 (2009.95). The most likely configuration is 
highly asymmetric, with a = 20.5^-3+o.s ^ A = Q> 70+0.03+0.05 ? 

and £ = -19°^ *2°> corresponding to a major-minor axis ra- 
tio of more than 2.4+q g, with symmetric models highly dis- 
favored. The resulting FWHMs of the minor and major axes 
are then 37 ± 1 //as and 88 ± 9/ias, though these are signifi- 
cantly correlated due to the substantially larger fractional error 
on A in comparison to that on a. The intrinsic image, scatter- 
broadened image and associated visibilities of the most likely 
configuration is shown in the middle row of Figure 3. 

The x 2 f° r m e combined data set is 61.5, and much lower 
than that found for the symmetric case. As described in Sec- 
tion 4, a decrease in x 2 is expected given the two additional 
addition of two parameters. However, the various ICs, given in 
Table 1, provide a means for identifying significantly lower x 2 - 
The best fit asymmetric model has a BIC that is 6.8 lower than 
the best fit symmetric model, and an AIC that is 10.4 lower than 



the best fit symmetric model. These provide "strong" evidence 
against symmetric models for the image of Sgr A*, ruling these 
out at 2.6a (BIC) and 3.2a (AIC) levels, in terms of the relative 
significance. That is, despite the limited visibility coverage in 
the u-v plane, the existing mm- VLBI observations can conclu- 
sively detect asymmetric structure in Sgr A*. 

5.4. Accretion Flow and Implications of Physics 

The images of the BL06 model, described in Section 3.1, 
are characterized by asymmetric crescents. These arise due 
to the combination of gravitational lensing, the relativistic or- 
bital motion and the opacity of the underlying accretion flow. 
The size and extent of the crescent depends upon the spin and 
inclination of the system, with nearly face-on disks (small 6) 
appearing annular. 

Because we have only have access to the visibility magni- 
tudes, configurations rotated by 180° are indistinguishable, im- 
posing an unavoidable ambiguity upon any results. In addition, 
despite opacity, the images exhibit nearly exact symmetry be- 
tween configurations viewed from above the equatorial plane 
(i.e., < 90°) and those viewed from below the equatorial 
plane (i.e., 6 > 90°) at equal inclinations. As a consequence, 
there is also an ambiguity in the line-of- sight component of the 
spin vector. For this reason, here we discuss only 0° < 6 < 90°. 
However we have performed the analysis for 90° < 6 < 180° 
as well, finding no statistically significant differences in the pa- 
rameter estimates. 

Figure 6 shows the minimum x 2 (maximum likelihood) as a 
function of spin, a, and viewing angle, 6. As with Figure 5, 
the four panels on the left show this for the individual epochs, 
while the large panel on the right shows x 2 f° r the combined 
data set. The morphology is similar in all cases, with the dif- 
ferent u-v coverage during the different epochs manifesting it- 
self primarily in the size of the likely region. As a result, we 
conclude that again the four epochs are consistent with a sin- 
gle underlying image morphology. Over all of the observations 
the flux normalization varied from 2.16 Jy (2009.96) to 3.00 Jy 
(2009.97), and is at all times sufficiently close to the value of 
2.5 Jy used to compute the images of the accretion flow. The in- 
trinsic image, scatter-broadened image and corresponding vis- 
ibilities of the most probable configuration (not necessarily the 
lowest x 2 , see Section 6) is shown in the bottom row of Figure 
3. 

At 53.09, the x 2 f° r the BL06 model (corresponding to a 
reduced-% 2 °f 0.830) is the smallest of any model we consider. 
Again we may assess the significance of the this by appealing 
to the ICs described in Section 4. In this case the BIC and 
AIC are given by 82.9 and 68.9, respectively. These are much 
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Figure 5. Asymmetric gaussian x 2 's as a function of a and A, at the most likely £ and Voo, for the 2007, 2009.95, 2009.96, and 2009.97 epochs, as well as when all 
epochs are combined. The minimum x 2 , number of visibility observations (N) and fit parameters (k) are listed in each plot. In all cases the color map ranges from 
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Figure 6. BL06 accretion flow x 2 's as a function of a and 6, at the most likely £ and Voo, for the 2007, 2009.95, 2009.96, and 2009.97 epochs, as well as when all 
epochs are combined. The minimum x 2 > number of visibility observations (N) and fit parameters (k) are listed in each plot. In all cases the color map ranges from 
min(x 2 ) (blue) to min(x 2 ) + 6 (red). 



lower than those from the symmetric gaussian, providing "de- 
cisive" evidence against symmetric configurations. Both are 
also below those of the asymmetric gaussian by 8.4 (since the 
number of parameters is the same in both models), implying 
"strong" evidence for the physically motivated accretion flow 
model in contrast to the phenomenological asymmetric gaus- 
sian. This corresponds to a 2.9a confidence level in terms 
of the relative significance of the two models; i.e., the phys- 
ically motivated accretion models are more than 67 times as 
likely as the phenomenological asymmetric gaussian models, 
and more than 2 x 10 3 times as likely as symmetric gaussian 
models. However, for two reasons this actually understates the 
case. 

First, the reduced-x 2 of the BL06 models is significantly less 
than unity. Indeed, with 64 degrees of freedom we expect 



a x 2 lower than that obtained (53.09) only 17% of the time. 
This suggests that the errors on the mm- VLB I visibilities have 
been over-estimated by a moderate amount. On the other hand, 
we may measure the "true" errors by assuming that the BL06 
model gives a sufficiently close approximation to the true im- 
age flux by renormalizing them until the reduced-x 2 * s unity. 
This requires a roughly 10% reduction in the errors quoted in 
Section 2. This in turn alters the x 2 , BIC and AIC for the other 
models as well. These values are listed under the "Implied 
Errors" section of Table 1. The net effect is to increase the sig- 
nificance with which the asymmetric gaussian is ruled out in 
favor of the BL06 accretion model to 3.2a. 

Second, the accretion flow models should generally be pre- 
ferred on the basis of their motivation alone. That is, there is a 
strong prior in favor of physically motivated models by virtue 
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Figure 7. Posterior probability of a given spin vector, marginalized over the Voo, obtained using the 2007, 2009.95, 2009.96, and 2009.97 epochs. Each sub-panel 
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of their design and connection to an existing body of knowl- 
edge. Furthermore, in the case of Sgr A* this includes the fact 
that these models were constrained to fit the pre-existing spec- 
tral data as well, something that is only possible because the 
physics governing the accretion flow provides a means to re- 
late the properties of images at different wavelengths. 

The fact that the BL06 model provides a significantly better 
fit to the mm- VLB I data implies that we can now distinguish 
phenomenological and physically motivated models on the ba- 
sis of mm- VLBI observations alone. Since the most prominent 
features of the accretion flow image are due to the generic prop- 
erties of black hole accretion flows, namely, the spacetime and 
orbital motion, this is likely to be robust among all radiatively 
inefficient accretion flow models for Sgr A*. That is, even with 
the extremely sparse u-v coverage presently available it is al- 
ready possible to probe signatures of general relativity and ac- 
cretion physics in the image itself. 

6. ESTIMATING BLACK HOLE SPIN 

Following Broderick et al. (2009) we produce posterior 
probability distributions for the parameters defining the BL06 



model, assuming a flat prior on a and isotropic priors upon 
the spin direction, and marginalizing over the Voo- For the 
3 -dimensional parameter space defining the vector black hole 
spin, (<2,#,0> this is done for each epoch individually (Figure 
7) as well as for the combined data set (Figure 8). For each 
epoch we show the probability distribution for the same posi- 
tion angle slices, defined such that the most probable values 
from the combined data set is exhibited in the central panel. In 
all cases we define la, 2a and 3a contours in terms of the cu- 
mulative probability, and these are shown by the solid, dashed 
and dotted lines, respectively. The probabilities are normalized 
to their average value, i.e., if all points were equally probable 
the probability density would be unity. 

Note that the way in which we have chosen the slices in £ in 
Figure 7 does not capture the most likely configuration based 
upon the 2007 epoch alone (shown in Broderick et al. 2009). 
Nevertheless, the regions shown are well within the la region 
from that epoch. Here we explicitly see that all epochs pro- 
duce consistent estimates for the spin, with varying degrees of 
statistical strength. 

The combined data set dramatically restricts the parame- 




Figure 8. Posterior probability of a given spin vector, marginalized over the Voo when all epochs are combined. Each sub-panel shows the posterior probability as 
a function of a and 6 at fixed £, set to the value shown in the lower-left corner of each. In all cases the center panel shows the most probable position angle when 
all epochs are combined, and the color scheme spans identical ranges. Contours demark the la (solid), 2a (dashed), and 3a (dotted) regions, defined by cumulative 
probability. 



ter estimates to a narrow sliver in the 3 -dimensional spin pa- 
rameter space. The most probable values for the spin are 

a = o.O +0 - 64+0 - 86 , 9 = 68°^ + ?; 8 o, £ = -52°%%, where the 
errors quoted are the la and 2a errors. At this point the proba- 
bility density is roughly 350 times that of the average value. In 
practice, these quantities are much more tightly correlated with 

<9 ~ 68° -42°a± 3° ±5° (15) 

Note that these are degenerate with solutions for which £ differs 
by 180°, i.e., £ = 128°!^o!^o, and for which the line-of-sight 

component of the spin is reversed, i.e., 6 = 1 12°_ 5 o_ 9 o . 

We do not attempt to determine the systematic uncertainties 
associated with selecting a particular accretion model. How- 
ever, we note that a number of efforts to fit alternative accre- 
tion flow models to the 2007 epoch have reached consistent 
results despite differences in the models, suggesting that these 
results are robust. Furthermore, the quality of the fits to the 
mm- VLB I visibility and spectral data, concurrently, suggests 
that the features of the BL06 model responsible for determin- 
ing the spectral and image properties are generic, and are there- 



fore insensitive to the particulars of the accretion flow models. 
However, full studies of the systematic errors associated with 
the particular choices made for the accretion flow properties 
and the underlying spacetimes are now justified. While beyond 
the purview of the present paper, we will report upon such ef- 
forts elsewhere. 

Probability distributions for each of the spin parameters, 
marginalized over all others, are shown in Figure 9. In addi- 
tion to the present case, these are also shown for the analysis of 
the 2007 epoch for comparison. In these the la and 2a ranges, 
defined by the cumulative probability, are also shown. In all 
cases the marginalized probability distributions from the com- 
bined data set are much more narrowly peaked than their 2007 
epoch counterparts. Nevertheless, they are all consistent at the 
la level with those obtained from the 2007 epoch alone. 

It is now possible to exclude a > 0.62 at the 2a level, with a = 
q+o.32+0.62^ substantially preferring non-spinning models. Thus 
the high-spin island seen in Figure 7 of Broderick et al. (2009) 
is now eliminated. Similarly, the position angle is now very 
clearly constrained, choosing the solution less favored by the 
2007 epoch data (though still consistent at the la level). In 
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Figure 9. Posterior probability of a given a (left), # (middle), and £ (right), marginalized over all other parameters, obtained using all epochs. The filled regions 
show the lcr (dark) and 2a (light) regions, defined by cumulative probability. 



this case we have £ = -60 o _ 5o _ 12 o . Finally, the most probable 



viewing angle is = 61 



o+7°+24 c 
-9°-15° 



. This is somewhat higher than 



the most likely value from the 2007 epoch alone, though well 
within the la uncertainty. 

These estimates for the orientation of the spin vector are in 
good agreement with a number of other efforts to estimate the 
properties of Sgr A*'s accretion flow. Estimates based upon 
fitting longer wavelength observations with numerical models 
of radiatively inefficient accretion flows produce position an- 
gles and inclination estimates with large uncertainties, though 
these are nevertheless consistent with the results obtained here 
(Huang et al. 2007). It is also in excellent agreement with more 
recent attempts to probe the spin orientation using the mm- 
VLBI data from the 2007 epoch (Huang et al. 2009; Dexter 
et al. 2010). As before it is not possible to assess consistency 
with models that employ qualitatively different plasma distri- 
butions near the black hole (e.g., Markoff et al. 2007), though 
they tend to imply similarly large viewing angles. 

We find similar spin orientations to those inferred from mod- 
elling of infrared polarization observations of Sgr A*'s flar- 
ing emission, though in this case the uncertainties are consid- 
erable (e.g., Meyer et al. 2007). Similarly, we find consis- 
tency with the spin directions obtained from modeling the spec- 
trum and polarization using general relativistic MHD simula- 
tions, despite preferring significantly smaller spin magnitudes 
(Shcherbakov et al. 2010). 

Unlike the estimates in Broderick et al. (2009), there is no 
longer any allowed solution for the spin vector that aligns with 
either of the reported stellar disks in the inner 0.2 pc of the 
Galactic center (Genzel et al. 2003). However, our revised po- 
sition angle estimates are consistent with being aligned with 
the X-ray feature reported in Muno et al. (2008), bolstering the 
interpretation of this as related to a possible jet. Note, how- 
ever, this interpretation may be inconsistent with the low spin 
magnitudes we prefer. 

7. OPTIMIZING FUTURE OBSERVATIONS 

The constraints upon the accretion model parameters ob- 
tained in the previous sections have implications for future 
mm- VLBI experiments. With these it is possible both to make 
predictions for the expected visibilities on the various possible 
baselines, as well as identify which baselines are most likely to 
provide substantial improvements to the BL06 model parame- 



ter estimation. To estimate these, here we compute the average 
visibility amplitudes as well as the variance associated with the 
uncertainty in the model parameters, weighted by the posterior 
probability distributions we have obtained using the combined 
mm- VLBI data set, following the method of Fish et al. (2009). 

The probability-weighted mean visibility profile of the 
scatter-broadened 230 GHz emission from Sgr A* is elongated 
in the (+w,+v) direction (Figure 10). For the moderately high 
values of favored by the mm- VLB I data, the intensity profile 
is dominated by Doppler-boosted emission on the approaching 
side of the accretion flow (e.g., in the northeast of lower left 
panel of Figure 3). This portion of the emission is elongated 
parallel to the projected direction of the black hole spin vec- 
tor. Since the effective size of the emission is larger along the 
projected spin axis, the correlated flux density falls off faster 
with baseline length for baselines that are sensitive to struc- 
ture in this direction than in the perpendicular direction. Our 
estimates of RIAF parameters suggest that a mm- VLBI base- 
line oriented southwest-northeast will detect more correlated 
flux density than an equal-length baseline oriented southeast- 
northwest. 

The standard deviation of the visibility amplitudes provides 
an estimate of which baselines would provide maximal addi- 
tional constraints on RIAF model parameters (assuming equal 
sensitivity at all sites). Previous computations based on the 
2007 epoch of data indicated that the largest scatter occurred at 
baseline lengths of approximately 3 GA and with orientations 
perpendicular to the Hawaii-SMT baseline (Fish et al. 2009). 
These findings still hold in light of the 2009 data. 

Among possible observing baselines in the next few years, 
Chile-SPT and Chile-LMT probe the region of highest standard 
deviation, followed by baselines between the LMT and the con- 
tinental US. Our model suggests that the LMT-SMT and LMT- 
CARMA baselines will detect well over 1 Jy. The increased 
sensitivity provided by phased ALMA may be important on the 
baselines to Chile, as the probability-weighted mean visibility 
amplitudes are > 0.1 /y on the Chile-SMT baseline, < 0.1 Jy 
over most of the (w,v) track of the Chile-CARMA baseline, 
and smaller still on the longer baselines to Chile. However, the 
standard deviation of the predicted model visibility amplitudes 
is several x lOmJy on these baselines, leading to uncertainties 
of tens of percent in model amplitude predictions. Further mm- 
VLBI data, either in the form of higher sensitivity on existing 
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Figure 10. Left: Mean predicted visibility amplitudes in Jy. Right: Standard deviation of predicted visibility amplitudes. In both cases the model image was 
weighted by its posterior probability. 



baselines or detections on new baselines, will both reduce these 
uncertainties and test the RIAF model with increasing rigor. 

8. CONCLUSIONS 

The significantly increased number of long-baseline visibil- 
ities, significantly larger signal-to-noise, and improved north- 
south coverage of the 2009 mm-VLBI observations have al- 
ready paid substantial dividends in the estimation of the prop- 
erties of Sgr A*. This is despite the fact that only three inde- 
pendent mm-VLBI stations (JCMT,CARMA,SMT) were used, 
and therefore the u-v plane remains extremely sparsely popu- 
lated, with long baselines primarily in the east-west direction. 
Constraints upon the black hole spin have improved dramat- 
ically in all cases, with the spin and viewing angle becom- 
ing tightly correlated, with the most probable configuration, 
a = 0.0 +0 - 64+0 - 86 , 9 = 6S°%t 9 _l S o , £ = -52°%% being roughly 
350 times as likely as the average probability density, and 25 
times as likely as the most probable configuration reported in 
Broderick et al. (2009). 

Despite the limited north-south coverage, the 2009 mm- 
VLBI data conclusively excludes symmetric gaussian models 
for the source. The preferred asymmetric gaussian has a major- 
minor axis ratio of 2.4^q ^q 3 6 , oriented with the major axis ori- 



ented -19 



o+3°+6° 
-l°-2° 



east of north. This implies major and minor 



axis FWHMs of 37 ± 1 /ias and 88 ± 9 //as, with the symmetric 
case excluded at 3.9a significance. Note that this orientation is 
not aligned with any particular feature in Sgr A*'s vicinity or 
the properties of the intervening interstellar electron scattering 
screen. 

It is natural to give physically motivated models a prior 
bias over phenomenological models. Nevertheless, even when 
physically motivated accretion models are compared with 
asymmetric gaussian models are weighted equally, the accre- 
tion models provide a significantly better fits to the mm-VLBI 
data. Based upon both the BIC and AIC we find strong evi- 
dence in favor of the accretion model, corresponding to a pos- 
terior probability 67 to 160 times larger than that of the most 
likely asymmetric gaussian model. This is particularly strik- 



ing given the simplicity of the accretion model we consider, 
suggesting that the image depends primarily upon the gross dy- 
namical and geometric properties of the system: the orbital mo- 
tion of the accreting material and the strong gravitational lens- 
ing by the black hole. In any case, it is clear that we have now 
entered the era of studying accretion and black hole physics 
with mm-VLBI. 
Our best fit accretion model requires a black hole spin of 



o+5°+9° 
-20° -28° 



oriented 



a = o.O +0 - 64+0 - 86 viewed at an angle of = 68 

at a position angle of £ = -52°*}5 *24°- Ambiguities due to 
the fact that only amplitudes of the visibilities were measured 
and the near symmetry of the accretion model images produce 
degeneracies corresponding to 6 <-> 180° - 6 and £ o £+180°, 
independently. The detection of closure phases will eliminate 
the ambiguity in position angle. 

During all of the mm-VLBI epochs presently available, 
three or fewer VLB I stations were employed, providing at 
most three baselines. In practice the two long baselines are 
nearly collinear, aligned predominantly east- west. The result- 
ing sparse coverage within u-v plane is the primary factor lim- 
iting the estimation of black hole and accretion flow parame- 
ters. Therefore, despite the success attained thus far, there is a 
considerable opportunity to improve the constraints upon Sgr 
A* and its accretion flow dramatically in the near future by 
including additional VLBI stations. The Event Horizon Tele- 
scope is a mm and sub-mm wavelength VLBI network whose 
goal is to observe, image, and time resolve structures near the 
black hole event horizon (Doeleman et al. 2009a). Over the 
next few years, new sites will join the current array, and sensi- 
tivities of critical baselines (e.g., those found in Section 7) will 
be enhanced through technical developments to widen band- 
widths and phase together multiple dishes at mm array sites 
(e.g., CARMA, SMA, ALMA). An important result is that 
within the next few years, interferometric phase information 
on long VLBI baselines will become available through robust 
measurement of the closure phase (Doeleman et al. 2009b), 
which will provide important new constraints on models of Sgr 
A*. 
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